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It is well known that uncorrelated, transverse displacements of the lenses 
of a beam-waveguide cause the light beam to deviate from its axis and that 
the tolerance requirements on the accuracy of the transverse lens positions 
are very stringent. 

This paper extends the statistics of beam waveguides to include correla- 
tions between displacements of different lenses and studies the effects of a 
succession of uncorrelated bends. 

It can be concluded from this work that the rms deviation of the light beam 
is proportional to the square root of the length of the waveguide if the corre- 
lation between lens displacements extends only over a limited range. 

The amplitude of the Fourier component of the waveguide axis whose 
period equals the oscillation period of the ray has to be less than 0.2 micron 
if the deviation of a light beam passing through 10,000 lenses of a confocal 
waveguide is to be kept less than 2 mm. This requirement means that the 
average radius of curvature of a waveguide composed of independent circular 
sections of an average length of 20 m with lenses spaced 1 m apart has to be 
more than 10 km. The comparison between two model guides, one composed 
of circular section and the other of sections shaped like sin 2 fix, indicates 
that the beam deflection depends only on the average radius of curvature and 
average length of the sections but not on their particular shape. 

I. INTRODUCTION 

Hirano, Fukatsu and Rowe 1 have studied the behavior of a light beam 
in a beam- waveguide whose lenses are randomly displaced from a per- 
fectly straight line. The first two authors considered also a waveguide 
with sinusoidal axis displacements. The behavior of light beams in bent 
lens-waveguides was studied in Ref. 2. These two papers represent two 
extreme cases of completely uncorrelated departures of the waveguide 
axis from perfect straightness on the one hand and perfectly correlated 
departures from a straight line on the other hand. 
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This paper describes the statistics of a light ray by introducing a cor- 
relation function connecting different points on the waveguide axis. 7 

We show how the ray position at the nth lens depends on one Fourier 
component of the curvature function of the waveguide axis, while the 
rms value of the beam displacement depends on one frequency com- 
ponent of the "power spectrum" of the curvature function. 

The dependence of the ray position on the Fourier component of the 
curvature function of the waveguide axis is analogous to the mode con- 
version loss of multimode waveguides. 3 

The description of the beam deflection in terms of the correlation 
function of the guide axis is used to draw some general conclusions. It is 
found that the rms value of the light beam deflection depends on the 
square root of the number of lenses in the guide provided that the cor- 
relation length is much less than the length of the guide. This fact can 
be used to deduce that the contributions of different, uncorrelated sec- 
tions of the waveguide to the mean square of the beam deflection simply 
add up, so that this value can be computed by computing the average 
value of the beam deflection of one section only. 

It is pointed out that apparently plausible models for the correlation 
function can lead to widely varying results. For this reason no attempt 
was made to describe the statistics of the waveguide in terms of the 
correlation function model. 

The results of this paper can be applied to alternating gradient focus- 
ing systems since it is known 4 that the beam deviation in such systems 
is of the same order of magnitude as that of a system of positive lenses. 

The results of this paper are of particular significance for beam wave- 
guides composed of gas lenses since such a waveguide would use closely 
spaced lenses so that the number of lenses for a given length of waveguide 
would be very large. The tolerance requirements are proportional to the 
square root of the number of lenses in the guide so that the tolerances of 
waveguides with closely spaced gas lenses become more stringent than 
those of waveguides using lenses spaced far apart. 

II. RELATION TO FOURIER SERIES 

We use the ray description of Ref . 2 to study the statistical behavior 
of a light ray in a lens-waveguide. 

The position of the light ray is given by its distance r„ from the lens 
centers. The inhomogeneous difference equation 2 

r n+ z - (2 — n)r n +i + r„ = Y n+2 (1) 

with 
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« = L/f (2) 

(L = lens spacing, / = focal length ) connects the ray positions at three 
successive lenses. The quantity F„+2 at the right hand side of (1) is the 
distance from the center of the (n + 2)th lens to the point at which the 
straight line through the centers of the nth and (n + l)th lens intersects 
the (n + 2)th lens (Fig. 1 ). If the lens spacing L becomes infinitesimal, 
L — ■> 0, Y/L assumes the meaning of the inverse radius of curvature R 
of the waveguide axis. However, even for finite lens spacing one can 
define a radius of curvature R n+l by the relation 2 (Fig. 1) 

L 2 

so that Y n is a measure of the curvature of the waveguide axis. 

The waveguide axis can also be described by the distance S n of the 
nth lens from an arbitrary straight line. 1 Between S n and Y n exists the 
following approximate relationship which can easily be derived from 
Fig. 1 

Y n +2 = — *5„ +2 "T" 2<S„+i — On . (4) 

The solution of (1) can be given in the form 2 

r n = r,/' + rj (5) 

with the solution of the homogeneous equation 

k „ , r\ — r cos 8 . ( . v 

r„ = j'o cos 11B H : — sin vB (5a; 

sin 8 

and the definition 

cos^ = 1 - \k (6) 

and the solution of the inhomogeneous equation 



H-i 
i 



-4- fl E Y«* sin (« - *)*, » ^ 2 - (7) 

sin 8 *=i 



We will calculate the ray's departure r„ from the axis at the end of 
the waveguide assuming that the ray entered the waveguide on-axis, 
r = n = 0. 

Equation (7) can he rewritten 

r„ = -7— Ism ti6 T] Y r+l cos v8 — cos nB >J F„ + i sin »>0> . (8) 
sin ff f=i *=i 
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Fig. 1 — Definition of various parameters of the beam -waveguide. 

The sums in (8) are Fourier coefficients as can be shown in the following 
way. The iV-discrete points Y y+1 can be represented by a Fourier series 

Kh-i = E «, exp [i(2w/NM V =1,2,-,N. (9) 

,i=-(JV-l)/2 

Equation (9) is a system of N simultaneous equations for N unknown 
quantities a„ . The Fourier coefficients a„ can be calculated by multi- 
plying (9) by exp [— %(2t /N)trv] and summing over virom 1 to N 



a, m I E Y r+1 exp \-i{2*/N)w]. 

IS y=l 



(10) 



The representation (9) works if N is an odd integer. Separating a, into 
its real and imaginary part 



a = a r + ?|8 S 



(ID 
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we obtain 

1 K 

a « = t? £ Y^+i cos (2tt/N)<tv (12a) 

0* = -iff E U sin (2ir/AV. (12b) 

iv *-i 

If we choose a such = 2irv/N (this may be possible only approximately 
but 2ira/N will approximate 8 closely if N is large) and denote the cor- 
responding values of a and by ae and /3j , we get from (8) 

n — 1 

r„ = — — -- S«s sin n0 + /3( cos rid}. (13) 

sin 

This is an oscillatory function with an amplitude r max which is a slowly 
varying function of n, 

W = ^^ V«« 2 + ft" = ^^ I a* I- (H) 

sin sin 

The maximum deviation of the ray in the vicinity of the nth lens is 
determined by the magnitude of the Fourier coefficient of K„ +1 belonging 
to the frequency 8. The amplitude \a§\ has to be extremely small to 
keep r mttX within reasonable limits. 
If S„ is strictly a sinusoidal function 

S, = A sin Bv (15) 

we obtain with the help of (4) 

Y, = kA sin (f - 1) (16) 

and from (10) with AT » 1 

| a, | = }$ kA. 

Equation (14) leads to 

nA(n — 1) 



2 sin 

Using (3) we can also write 

L\n - 1) 
2 B sin 



(17) 



' nifkx — 



where R is the minimum radius of curvature. 

If we take k = 2, L = 1 m, n = 10,000 we obtain a waveguide of 10 
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km length. Requiring r max % 2 mm we find that the amplitude A, the 
maximum departure of the waveguide from a straight line, has to be 

A ^ 0.2 micron 

or 

R = 2500 km. 

These numbers show how extremely small the 0-Fourier component has 
to be! 

If we take Y,+i = B sin &v and substitute into (10) with 8 = (2r/N)a 
we get nonvanishing values for | a B | even if # ^ 8 and TV — > so. These 
nonvanishing values appear for all # values satisfying the equation 

j? = 8 + 2*? (p = 0,1,2,3, - - • ). 

It appears, therefore, as if we obtain Fourier components a s for all 
harmonics & = 6 + 2irp. This apparent discrepancy is resolved if we 
consider the period length X,? = 2irL/& of the oscillation sin ($/L)vL. 
The above equation leads to the solutions for the period length 

L\b fa if p = 

Xtf = 



L + pXa \<L if V ^ ° 

with X = 2rL/8. The period length Xd is therefore either equal to X e , 
the natural period of the ray oscillations, or it is less than h. Since the 
lenses are spaced a distance L apart a period of length less than L is 
meaningless! 

III. RANDOM DISPLACEMENTS OF THE GUIDE AXIS 

Our considerations so far have been limited to a definite shape of the 
waveguide axis. However, they can easily be extended to a statistical 
theory. Equation (14) can be used to obtain the rms value of the maxi- 
mum beam displacement. 

a = v^ = ^v<rwF). (is) 

The symbol { ) designates an ensemble average. The quantity ( | a B | 2 ) 
is the expected value of the 6 component of the "power spectrum" of 
the waveguide curvature. It is also possible to express A not by means 
of the power spectrum but by the correlation functions of Y, . For this 
purpose we write, with the help of (10), 

(I a, f) = 4 Z £ (Y„ +1 Y^) exp [0(ji - p)] . (19) 
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It is reasonable to assume that (F p +iF„ + i) depends only on the difference 
M — v so that we can write 

<F, +l 7 M+1 > = /„_, (20) 

and to assume, furthermore, 

Jfi— v = J v— M ■ ("1/ 

The factor /x is the correlation function of the curvature function of 
the waveguide axis. 
Equation (19) can be rewritten in the following way: 
, r i(jv-i) ar — 1*| 

(|a.| ! ) = p £ E Ue™ 

N* {, l(N-l) t-\*\+l 

+ E E /2.-1 exp [#(2s - 1)] 
#-1 t-i«i+i 

-l(AT-l) W+l-]i| 

+ E E /an-i exp [#(2s + 1)] 

The summation over t can be carried out since the terms under the sum- 
mation signs are independent of /. Using (21) to simplify our expression 
further we obtain 

1 ! (Jv-i) ^ 

(22) 



so 


that 


(1 «• I 2 ) = 
(18) becomes 


1 c (Jv-i) 

*{** + »£ 


(JV - 
X 

n — 


X)/x 


cos 

COG 


xal 




A- Vn ~ 

sin 


iX0>. 



(23) 

If /x decreases with increasing X so that the upper limit in the sum of 
(23) becomes immaterial the equation shows that A is proportional to 
Vn-1 and not to n — 1 itself as one might have suspected by looking 
at (18). If the waveguide's curvature contains a sinusoidal component 
which persists throughout its length, so that no correlation length less 
than n exists, then the sum is proportional to its upper limit n and A 
becomes truely proportional to n itself. 1 Equation (23) expresses the 
rms-beam deviation in terms of the correlation function /x of the curva- 
ture of the waveguide axis. It is easy to rewrite this equation as an ex- 
pression depending on the correlation function of the waveguide axis 
displacement itself. The displacement of the nth lens from a straight 
line is S n . Defining 

{S V S») = <?,_, = (?„_, (24) 
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we get with the help of (4) and (20), 

/ x = 6(r\ - 4G X -i - 4CX+1 + Gx-2 + &+s . (25) 

Substituting (25) into (23), rearranging terms leads to 

A - Vn- i« | Go + 2 f[i_ X/(n- DlffxoosXfl 

sin 6 ^ x=i 

+ 7 1— [2(1 + 1 - fctyft - (2 - c)ft 

(n - IV L 

00 

- 4k V* - i* 2 E ^x sin X0 

x=i 

where we assumed that the correlation length A « n. If we assume that 
n » 1 (26) simplifies 

A a V5« (g 4- 2 £ Cx cos XflV- (27) 

sin [ x=i J 

This approximate equation shows that A is very nearly proportional to 
the 8 component of the Fourier coefficient of the correlation function. 
If the lens displacements S, are un correlated, Gx = for X ^ 0, we 
have 

A = 2 l/jJL-aVft (28) 

with 

8 = VGo = VW>* 

Equation (28) differs from Ref . 1 (17) (for large values of n) by a factor 
of \/2- The reason for the occurrence of this additional factor in our 
theory is that our A is the rms value for the amplitude of the oscillatory 
beam trajectory while (17), Ref. 1 describes the rms value of all points 
of the oscillatory beam trajectory. 

Let us assume that the waveguide axis is composed of sections of a 
given average length and that the curvature function of one section is 
uncorrelated to that of any of the other sections. All sections are assumed 
to be of the same type. For example, they may all be circular bends which 
differ only in length and radius of curvature. A waveguide of this type 
has a finite correlation length and if n, the total number of lenses, is 
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large we get from (23 ) 

A n 2 = An (30) 

with 



A = 



(si^y^ +2 s^ c ° sx i (3i) 

Adding one more section with m lenses changes n to n + m in (30). The 
increase of A 2 due to the addition of a section with m lenses is given by 
A m 2 = Am so that we get 

A n+m 2 = A n 2 + Aj (32) 

for a guide with n + m lenses. 

Each section of the guide can be thought of as such an additional sec- 
tion so that 



^ = Z < 2 



s 

i=l 

with A„ r 2 being the contribution of the *>th section and M the number of 
sections. Introducing the average value A m 2 , wo obtain 

A = A m VI. (33) 

The rma beam deviation of a waveguide composed of sections can be 
obtained by calculating the contribution of the nns beam deviation of 
each section, computing their nns values and applying (33). 

IV. EXAMPLES 

As a first example we consider a waveguide composed of circular sec- 
tions which are connected so that the first derivative is continuous. The 
departure of the light beam which goes through the center of the first 
two lenses of one of the circular arcs is given by (7) 

sin n = sin (n — 1) - 
r n = Y - (34) 

a 

sin sin - 

with Y = Y „ = const. If the circular sections have an average number 
of m lenses which vary with a Gaussian distribution with variance <r m we 
get with the help of (33) 
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A = vmmL +icose 



sin 8 sin - 



(35) 



2 cos - cos (m — \) 8 exp (|0 a m ) 



+ | cos (2m - 1)8 exp (-20*0 



with M circular sections per waveguide. It was assumed that Y and m 
are statistically independent. 

The validity of (35) was checked by a computer simulated experi- 
ment. We constructed 30 different waveguides composed of a series of 
circular arcs. The quantities Y and m were computed as Gaussian ran- 
dom variables with mean value { Y ) = and ( m ) = m. Each wave- 
guide contained 10,000 lenses. One light ray was traced through each 
guide with the use of (1 ) and the rms value of the values r n with n = 
10,000 was computed which, multiplied by \/2, should equal the value 
A of (35). This experiment was repeated three times with different 
values of m. In all experiments we considered the confocal case, k = 2, 
cos 8 = 0. 

Table I shows how (35) compares to the computer results. The a m 
values in Table I were chosen for the convenience of the computer 
calculations. 

The agreement between the theoretical and experimental values is 
quite good considering that the rms value was computed from only 30 
samples. 

We can relate { Y 2 ) to an average radius of curvature R since ac- 
cording to (3) 



<n = l 



1 \r 2 / r 2 



(36) 



Table I 



» = 10,000 



m 


"m 


A 
, (equation (35)) 


Computer 

y/T^i) Experiment 


3 

20 

100 


0.92 
6.1 
30.0 


77.9 
22.4 
10.0 


00.2 

21.9 
8.7 
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Table II 
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m 


R 


3 

20 

100 


39 km 

11.2 kin 

5.0 km 



We may allow A = 0.2 cm at the end of the waveguide of 10 km length 
(L = hn). The pcrmissable values of R, computed from the theoretical 
values of Tabic I, are listed in Tabic II. 

As a second example we consider a guide which is built up of sections 
of tapered bends formed according to Fig. 2 



S, = A sin 2 0r/?n)O - 1) v = 1,2 



m. 



(37) 



The amplitudes A and the number of lenses m are random variables. 
Substituting (37) into (4) and (7) we obtain 



?*„+i = — A 



■ 2 T 
Sill — 

m 
sin 8 



(38) 



( sin On ) ( cos cos 6 j + ( sin d\( cos — n — cos Bn j 

2 sin (I " l) s!n (^ + 1) 

The lens numbered n = m + 1 is the last lens on the bend. According 
to (5a) the amplitude of the oscillations is given by 



= r " + 



( n — Tp COS g \" 
sin 6 )' 



(39) 



We restrict ourselves to the confocal case k ~ 2, cos 8 = and set r Q = 
fm , T\ = r m+i . because the oscillation caused by the bend is taken as the 




Fig. 2 — Beam -waveguide composed of sections of tapered bends. 
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initial condition of the ray. We obtain from (38) and (39) 

(sm^Vl + cos^Vl - cos^m) 
rm J = A * \ ™/\ EZA L-l. (40) 

■["(i-y'^s+i)*] 

To compute A we must average r max 2 over A 2 and m. We assume that 
A and m are uncorrected and that m follows a Gaussian distribution 
with average value m and variance a m . To simplify the averaging trim 
is replaced by w/m and these terms are taken out of the integral. This 
procedure is valid only if irlm « 1 so that we may also take sin ir/m = 
ir/m and cos 27i7m = 1. Remembering that M = n/m, we obtain from 
(33) 



VW) 



ir 2 \/M A/ l — cos ^ m exp ( — o n 2 * 3 '™ ) 



(41) 



m -sm 



with n being the total number of lenses of the guide and m the average 
number of lenses per section. 

The equation (41) is valid if m » 1 and tr m /m « 1. It shows that 
A decreases as the number of lenses per bend is increased. 

Some numerical results are shown in Table III. The column labeled 
"computer result" again contains data calculated by tracing rays through 
simulated beam-waveguides. The waveguides were "constructed" of 
arcs according to (37) with A and m being Gaussian random variables. 
Thirty random waveguides were constructed for each value of ni and 
the rms values of the ray position r n with n = 10,000 was computed from 
these 30 samples. 

The computer results agree with the theoretical values to the order 
of magnitude. The agreement in this second example is poorer than that 
of Table I. However, the 30 values of rio.ooo of this example scatter more 
widely than those of the first example. Omitting the largest of the 30 
values for m = 50 changes A by a factor of 0.65 which shows that the 
statistics of these 30 samples is not very reliable. By comparison omitting 
the largest value of the sample with m = 20 of the first example changes 
A only by a factor of 0.95. 

The results of Table III can be used to get an impression of the line 
tolerances required for a nominally straight waveguide. Let us assume 
that the 10,000 lenses of our model guide are spaced 1 m apart (L = lm) 
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Table III 
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n = 10,000 



m 


<r m 


A 
,. equation (41) 


Computer Result 


10 

50 

250 


2.8 
15 
77 


7.70 
1.12 10-> 

2.00 10- J 


12.6 
8.8 10" 1 
2.7 10- 3 



which results in a guide 10 km long. If this is a guide built of gas lenses 
of 0.05 cm diameter (about 4 in.) we might allow A = 0.2 cm to be 
reasonably sure that the light beam will get through the pipe. Fixing A 
allows us to calculate the rms value of the amplitudes -\/( A 2 ) of the 
deviation from straightness and the average radius of curvature of the 
guide. An average radius of curvature R at the peak of the arc of (37) 
is given by 



1_ 
R 



m-L- 



21 



V: 



m 



■\/n L 



(i-sJ'^S + ii)'- 



(42) 



sin l - — — 



The second half of this equation was obtained by substituting (41 ) for 
(A 2 )* assuming that \-kgJ » 1. Table IV shows the values of the rms 
amplitudes and the average permissable radius of curvature which were 
computed from the theoretical values of Table III. The tolerance re- 
quirements arc rather stringent as Table IV shows. For short bends, 
departures from straightness of only a fraction of a millimeter can be 
allowed while this tolerance moves up into the meter range as the aver- 
age length of the bend exceeds hundreds of meters. In case of uncor- 
related random wiggles our present example of A = 2 mm and n = 
10,000 lenses leads to an rms value for the position tolerances of 5 = 
0.01 mm according to (28). This is a very real tolerance requirement 

Table IV 



m£ 


VU$ 


s 


10 in 

50 111 

250 m 


2.0 lO -8 cm 
1.78 em 
1.00 10* cm 


19.4 km 
7.1 km 

3.16 km 
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since there will always be random lens displacements superimposed on 
longer bends so that the uncorrelated, random component of lens 
displacements has to be kept below 10 microns. 

A comparison of Table II with Table IV shows that the average radii 
of curvature permissable for our two examples are nearly the same. This 
is true even though the waveguides of the two examples are of very 
different construction. This may support the belief that the average 
radius of curvature and the length of the sections determine the deflec- 
tion of the ray regardless of how the waveguide is shaped in detail. The 
agreement between the values of Table II and IV is improved if one 
corrects for the different length of the sections. 

It may be in order to add a remark concerning models for the correla- 
tion function G\ , Since the correlation has to be of finite length one 
might be tempted to try a correlation function of the form 

G* = (? exp(-i^J) (43) 

with q being the number of lenses within the correlation distance. The 
correlation number q must be of the same order of magnitude as the 
average number m of lenses per section of our model waveguide. Sub- 
stituting (43) into (27) we obtain 

Bin 6 __ (44) 

a A V ( -i / > ex p t- 1 /^ - f;ose 

-y 1 - 2 exp {-l/q) l + exp ( _ 2/g) _ 2 exp ( _ 1/ff y COB 

or if l/g«l 

A = Vnjc A A Go . (45) 

sin0 T q(l - cosfl) 

This correlation function is obviously a poor model for a waveguide 
with random bends since A of (45) decreases with q only like q~* while 
A of (41) decreases like m~ 2 . For q — 250 we obtain from (45), setting 
VGo & y/W) (G » hVW) for the arcs of < 37 ))> 

which is three orders of magnitude larger than corresponding values of 
Table III. 

Another possible choice for a correlation function may be 

Gx= <7 exp (-X'/fl 1 ). ( 46 ) 
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Using the identity 5 

QQ DO 

£ exp (-X7<T) cos \9 - v^ g £ ex P (-t¥[(0/2t) - »] 2 ) (47) 

we obtain from (27) 

G~^r v * f exp i-'Vi(^) - if) 

Lro sm p (48) 

+ exp (-i?V) + exp (-7rY[(*Ar) + if) |* 

where all but three terms of the sum on the right hand side of (47) are 
neglected which is justified if q > 1, 
The maximum rms beam deviation of (48) decreases with increasing q 

like 

V^exp (-£?¥) 

that is much faster than (45). These two examples indicate hSw criti- 
cally A depends on the shape of the correlation function. It appears that 
more insight can be gained by choosing models for the random deviation 
of the waveguide curvature rather than by trying to guess at model cor- 
relation functions. 

The reason for the critical dependence of A on the shape of the corre- 
lation function can be seen from the following argument. 

In our second example, (37) we obtain Go ^ 0.5 (A 2 ). However the 
ratio A/y/(A 2 ) is, for example, in the order of 0.1 according to the 
second line of Table III. From (27) we obtain, with n = 2 and n = 10 4 , 

If A/\/(A-} is to be 0.1, the sum under the square root sign must be 
very nearly equal to —0.25 so that the two terms under the square root 
sign cancel to a term of the order of magnitude 10 -6 . A very slight varia- 
tion of the value of the sum gives rise to a large variation of A. 

As a last example we consider a waveguide with M random tilts. If 
each tilt is located at a lens, Ref. 2 (31), gives for the beam amplitude 
caused by one tilt with angle a 

(49) 



V4k - 
so that (33 ) leads to 



A = 2^/(rf)—7==^VM. (50) 

V4k — /C 
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If 7i = 10,000, L = lm, k = 2, and M = 100 we find that 

V{c?) S 2 10^* radians = 0.0115° 
if A ^ 2 mm is required. 

v. CONCLUSION 

Is it more advantageous to space the lenses of a beam-waveguide 
closely or farther apart? The answer to this question depends on our 
ability to control tolerances. Equation (41) shows that the mis beam 
deviation due to random bends decreases rapidly with increasing num- 
ber of lenses, while (28) indicates that the rms beam deviation due to 
uncorrected lens displacements increases slowly with increasing lens 
number. Only practical experience can tell how to compromise between 
these two conflicting requirements. 

APPENDIX 

Equivalence of Two Representations 

The problem of ray propagation in a bent lens-waveguide has been 
treated in Ref. 2 and by Hirano and Fukatsu. 1 The treatments of the 
problem in these two papers differ in the way the ray is described. In 
Ref. 2 the ray position r„ is measured from the center of the lenses and 
the position of the lenses with respect to each other is described by a 
quantity Y n (Fig. 1). In Ref. 1, a straight reference line is used to 
determine the position p„ of the ray as well as the lens displacements S n . 
The ray position r n at the wth lens in the representation of Ref. 2 is 
given by 

1 n-l 

r n = ~r— n E *Vi sin (n - v)B n ^ 2. (51) 

sin ?=i 

The values of r„ at n = and n = 1 are r = n = 0. 

In the representation of Ref. 1 the solution of the inhomogeneous 
difference equation reads 

Pu = -As E S, sin (n - p)$ n^2 (52) 

sin d y=\ 

with po = pi = 0. 

With the help of Fig. 1 it is easy to see that approximately 
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S n+2 = 2S n+1 - S n - Y, i+ iVl - (S n+l - S n )*/L*, (53) 



or if 



z 

Y n+ 2 = 2>S'„ + i — S n+ 2 — S n , (53a) 

and 

Pn = r n + S n . (54) 



The substitution of (53a) and (54) into (51) leads to 

-i n — 1 

Pm = S n -\- - — - zl (2£„ — S v+ i — S„-i) sin (n — v)d 
sin 8 t={ 

= S„ + -r— n {Si sin (n - 1)6 - S n sin 9 - So sin (n - 1)8 
sin 8 { 

n—l 

+ S *S„[2 sin (n - i>)0 - sin (n - v + 1)9 - sin (« — v — 1)8] 

which can be simplified to 

p„ = (So — Si) cos nd — (So — Si) cos 8 sin n0 

+ -r^ E & sin (n - ?)ff. (55) 

sin r=f 

The first two terms with *S — & are solutions of the homogeneous dif- 
ference equation and appear here because r = ri = does not coincide 
with pn = pi = if So and &i are unequal to zero. 

Since the first two terms can always be removed by adding a suitable 
solution of the homogeneous equation the equivalence of (51) and (52) 
has been shown. 
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